No spin-localization phase transition in the spin-boson model without local field 
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We explore the spin-boson model in a special case, i.e., with zero local field. In contrast to 
previous studies, we find no possibility for quantum phase transition (QPT) happening between the 
localized and delocalized phases, and the behavior of the model can be fully characterized by the 
even or odd parity as well as the parity breaking, instead of the QPT, owned by the ground state 
of the system. Our analytical treatment about the eigensolution of the ground state of the model 
presents for the first time a rigorous proof of no-degeneracy for the ground state of the model, which 
is independent of the bath type, the degrees of freedom of the bath and the calculation precision. 
We argue that the QPT mentioned previously appears due to incorrect employment of the ground 
state of the model and/or unreasonable treatment of the infrared divergence existing in the spectral 
functions for Ohmic and sub-Ohmic dissipations. 
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A two-level system coupled to an environment pro- 
vides a unique test-bed for fundamental interests of quan- 
tum physics. Denoting the environment by a multi- 
mode harmonic oscillator, the spin-boson model (SBM) 
|l|,l2J presents a phcnomcnological description of the open 
quantum system, which plays an important role in quan- 
tum information science and condensed matter physics. 
Particularly, for the SBM at zero temperature, it has at- 
tracted intensive interests for the quantum phase transi- 
tion (QPT) happening between localized and delocalized 
phases regarding the spin. 

The standard SBM Hamiltonian in units of K = 1 is 
given by [l[, 
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ui k ala k +^X k (al + a k )a z , (1) 



where a z and a x are usual Pauli operators, e and A are, 
respectively, the local field (also called c-number bias [1() 
and tunneling regarding the two levels of the spin. a\ and 
a k are creation and annihilation operators of the bath 
modes with frequencies uj k , and X k is the coupling be- 
tween the spin and the bath modes. The effect of the har- 
monic oscillator environment is reflected by the spectral 
function J(w) = 7r^ fe ^^(^ — w fc) for < ui < ui c with 
the cutoff energy lo c . In the infrared limit, i.e., u> — ^0, the 
power laws regarding J(o>) are of particular importance. 
Considering the low-energy details of the spectrum, we 
have J(w) = 2irau>l~ s uj s with < lu < lu c and the dis- 
sipation strength a. The exponent s is responsible for 
different bath with super-Ohmic bath s >1, Ohmic bath 
s =1 and sub-Ohmic bath s <1. 
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The local field e introduces asymmetry in the model, 
which was considered to be less important than the 
tunneling A and thereby sometimes neglected for con- 
venience of treatments. For the Ohmic dissipation, it 
was mentioned 0, [j that the SBM has a delocalized 
and a localized zero temperature phase, separated by a 
Kosterlitz-Thouless (KT) transition in the case of e = 0. 
In the delocalized phase, realized at small dissipation 
strength a, the non-degenerate ground state behaves like 
a damped tunneling particle. In contrast, for large a, 
the dissipation leads to a localization of the particle in 
one of the two a z cigenstates, implying a doubly degen- 
erate ground state. On the other hand, intensive studies 
have recently been paid on the sub-Ohmic dissipation, 
which also demonstrated the QPT between localized and 
delocalized phases [3J-422I] . In particular, the QPT could 
also happen in the absence of the local field (e = 0) [7H9| 
from a non-degenerate ground state with zero magneti- 
zation below a critical coupling to a twofold-degenerate 
ground state with finite magnetization for a coupling 
larger than the critical value. The QPT has been found 
so far by different numerical studies, such as numerical 
renormalization group (NRG) technique 0, @-|a IfOrJLT 
quantum Monte Carlo [9[ , the exact diagonalization 
the density matrix renormalization group approach 
and variational matrix product state approach [f9|. 
alytical studies, such as unitary transformation method 
[2(1 |21| . were also employed. The latest work by varia- 
tional method also showed the appearance of the QPT 
between delocalization and localization |22| in the case 
of e = 0. 

Concentrating on the SBM without the local field, 
we argue in the present work no possibility with the 
QPT happening between the localized and delocalized 
phases in above dissipative cases, opposite to previous 
results. The main behavior of the model could be char- 
acterized by the parity keeping or breaking, rather than 



the QPT, owned by the ground state in the variation of 
some characteristic parameters. The remarkable feature 
of our treatment is the analytical investigation, which 
is independent of the bath type, the degrees of free- 
dom of the bath and the calculation precision. The key 
step of our treatment is the decomposition of the spe- 
cial SBM Hamiltonians in such case into two decoupled 
sub-Hamiltonians, and the key point is the focus on the 
parity of the ground state of the model. We show that 
the significant changes regarding the ground state of the 
model, if can be called 'QPT', happen only when the par- 
ity breaks down, corresponding to appearance of the local 
field. Our investigation of the magnetization analytically 
indicates that the QPT discussed previously by differ- 
ent numerical studies is due to unreasonable treatment 
of the ground state and/or of the infrared divergence in- 
herently existing in the spectral functions for Ohmic and 
sub-Ohmic dissipations. 

We get started from the special SBM with zero local 
field. In such a case, we denote the Hamiltonian H by 
Hsb with 
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SB 
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a k + ^2^k(al + a k )a z , (2) 
k 



and introduce a parity operator 



(3) 



to commute with Hsb, i.e., [Hsb,^1] = 0. Moreover, by 
introducing a unitary transformation 



U = 



1 




odd-parity with respect to the parity operator a z since 

Meanwhile, a z = UHU + is actually the parity opera- 
tor of Hsb after the unitary transformation U has been 
performed, which means that Hsb also consists of two 
sub-Hamiltonians with the eigenfunctions of even-parity 
and odd-parity, respectively, with respect to the parity 
operator II. 

If QPT happens in the model of Hsb, there should 
be the possibility of degenerate ground states of the two 
decoupled sub-Hamiltonians H + and H~ and also the 
possibility of releasing the degeneracy @, 13, 0t11- To 
check the possibilities, we may solve H + and H~ col- 
laboratively by expanding \(p + ) and \ip~) using displaced 
coherent states 
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where cf n} 



are coefficients to be determined, {n} = 

tin arc for different bosonic modes, ({ti}}^ = 

,t 



K + qk) nk e~ q " a k |0) is the product of dis- 



placed coherent states of different modes of the bosonic 
field with the displacement variables q k = X k /oj k , k = 
1, 2, • • • ,N, and |0) the vacuum state of the bosons. 

Expanding the tunneling terms in Eq. ([5]) by {|{ra})^} 
yields 



v- t A 
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E D {m,n} \{m}) AA {{n}\ , (7) 
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{m},{n} 

l} is given by [Hj 



we may diagonalize Hsb to be two decoupled sub- 
Hamiltonians, 



H d = UH SB U + = 



H+ 
H~ 
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where 



H+= -_ e <* £*«!»* + Y^Aa k + A fe (4 + a k )], 

k 
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(5) 



It implies that the eigensolution of Hd consists of two 
independent eigcnsolutions of H + and H~ . Since no lo- 
cal unitary transformation can change the eigenvalues, 
the ground state of Hsb should be the lower one of the 
ground states of H + and H~ . 
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For the eigcnsolutions Hd I ' _ I = E + I n 1 and 
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Hd ( i -\ I = E I , 1, I, it is easy to prove that 
the eigenfunctions are, respectively, of even-parity and 
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So the prefactor e~ 2 ^fe 9 fc, relevant to Di mn \ and the 
tunneling terms, plays a crucial role in the SBM and 
should be non-zero. Otherwise the SBM in our case is 
physically meaningless due to neither local field nor tun- 
neling involved. This point will be further emphasized 
later. 

The solutions of H + and H~ by Eq. (6) yield 
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Thus the lowest eigencnergics for the even- and odd- 
parity states are, respectively, 
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It is intuitive to reach degeneracy by setting A = 0, which 
is, however, physically meaningless in our case. For any 
non- vanishing value of A, we show a rigorous proof in 
Appendix that there is no possibility of £t, = E 



{<>} 
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That is to say, a physically meaningful solution of the 
SBM in our case requires no degeneracy for the ground 



states of 



¥> n 



and 



<P 



which implies no QPT happen- 



ing, in terms of the viewpoints in [iMa, LZH9|, no matter 
how to change the characteristic parameters and which 
bath is considered. 

Our results also mean that the ground state of Hsb is 
surely of a certain parity, and any employed state with 
mixed parity (e.g., the superposition of the ground states 
of \<p + ) and \(fi~)) is definitely not the ground state of the 
model. In this sense, it is not a surprise that QPT was 
found in [221 ] using a presumed ground state which, due 
to mixture of odd- and even-parity, is actually not for 
the ground state of the model. In fact, our result could 
be obtained more simply by calculating the magnetiza- 
tion M = (cr z ). Without losing generality, we assume 
the eigenfunction of Hd with the two sub-eigensolutions 

. J _( ), where lip^) are the 
,sm% )y' '^ 0/ 

ground states, respectively, of ly*). So the magnetiza- 
tion of the Hamiltonian Hsb is 



correlated, i.e., |\fr) 



M = {a z ) = (V\Ua z U + \V) = -sin2(%J|e" r ^ Q 'H^) 

(11) 
If 9 = k-K or (k + l/2)ir with k = 0, 1, • • • , which corre- 
sponds to even-parity case or odd-parity case, the mag- 
netization M is absolutely zero, which implies no QPT 
happening. The finite value of M occurs only in the case 
of the superposition of \pq) and \<Pq), which, according 
to our discussion above, is not the ground state of the 
model. 

Besides the incorrect employment of the state as the 
ground state [22| , the main reason for obtaining the QPT 
in previous literatures |3l— [22| is the unreasonable treat- 
ment of the infrared divergence existing in the spectral 
functions for Ohmic and sub-Ohmic dissipations, which 
induced the degeneracy of the ground states in the low 
frequency domain. To clarify this point, we demonstrate 
below the potential infrared divergence in the treatments 
using the bath mode with continuous and discretized 
spectra, respectively. The NRG is taken as an exam- 
ple to show that the prediction of the QPT in previous 
studies is completely due to divergent expansion in the 
infrared limit. 

For treatment with the bath mode of continuous spec- 
trum: We return to the prefactor e -2 ^ fc 9fc , for which, 
using the spectral function J(ui), we have, 



with 

^[1 -(£)*-] if S <l 

ft =< Ms?) ifs=1 

where a and w c are defined above in the spectral function. 
ui± is a small quantity regarding the frequency difference 
from u = 0. In the case of the infrared limit, i.e., Wi — > 0, 



have J2k1k ~* °° ^ s ^ 1> which leads to D 
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J {o} = E {o} from E * ( 10 )' 



(See Eq. (8)) and thereby £7 

As a result, with calculation approaching the low fre- 
quency domain, the NRG gives a significant change from 
non-degeneracy to degeneracy, which was called QPT. In 
addition, we see from Eq. (7) that the tunneling terms 
arc vanishing with the calculation approaching the in- 
frared limit, which reduces Hsb to be 



#' = £< 



u k a\a k + ^ ^*(°l + a k )a z , 



commuting with a z . This corresponds to the localiza- 
tion, sometimes also called the frozen spin [8|. It means 
that, once the system falls into a well-polarized state of 
the spin, the magnetization will remain unchanged in the 
dynamics. Consequently, it is the infrared divergence 
that wrongly predicts the QPT happening from the de- 
localization to the localization as the infrared limit is 
approached. 

Eq. (12) presents the relationship between the pref- 
actor and the dissipation strength a, from which we see 
no possibility to have a finite value of a as the criti- 
cal number for the QPT. The degeneracy of the ground 
states occurs only in the case of a — > oo, i.e., a physi- 
cally meaningless case. On the other hand, considering 
Eq. (7), we may find that the infrared divergence also 
yields the parity operator II = o- x e m ^* a k a *> to be zero. 
In fact, the correct understanding of the parity variation 
could be reached only after we have eliminated the in- 
fluence from the infrared divergence. In such a case, the 
parity of the ground state of the model will explicitly 
break down as long as we introduce the non-zero local 
field into Hsb, which brings about a significant change 
in the magnetization as a function of e. We prefer to call 
this phenomenon as parity breaking, instead of the QPT. 

For treatment with the bath mode of discretized spec- 
trum: By means of the NRG logarithmic discretization 
parameter A defined in [4|, 0], we have £ fe Ik — 2a/?2, 
where 
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with N a large number required by the method [j] , (3q = 
( s + 2) 2 (l-A-s- 1 ) 3 /[( S + l) 3 (l-A- 8 - 2 ) 2 ]. So £*«£-► 
oo in the case of N — > oo and s $J 1. As an example, 
we show in Fig. 1 the distribution of /3 2 with respect 
to the mode- number N for different s. The divergence 
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FIG. 1: (Color online) Divergence of P2 in the calculation by 
NRG, where A = 2, w c = 1, and s = 0.1 (s = 1.0 in inset). 



and would be helpful for our deeper understanding of the 
weird phenomena in open quantum systems. 
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Appendix: Proof of the impossibility of degeneracy 

The proof starts from Eq. (10). Excluding the trivial 



case A 



0, we have Et Q -, 



B fnl occurring only for 
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(!) J2{n} C tn} D {0,n} = Or (2) £ { „ } C+ n} £» {0i „ } /c+ 
-T,{„}C{n} D {0,n}/ C {0} ^°- 

Employing the result 



D {0 . n} =e- 2 ^^l 



(2q k y 



leads to E~Ly = -Er i in the case of Ohmic and sub- 
Ohmic dissipations, from which the QPT was predicted 
in previous studies. 

Our work not only indicates for the first time that the 
parity keeping and breaking, rather the QPT, can fully 
characterize the behavior in the SBM, but also helps to 
understand a fact that the QPT had only been found 
by previous numerics in the sub-Ohmic and the Ohmic 
cases, instead of in the super-Ohmic case, which is be- 
cause the infrared divergence never happens in the super- 
Ohmic case, e.g, ^2 k q\ = 2a/ (s— 1) for s > 1 in the case 
of the continuous bath-mode spectrum. In fact, our in- 
vestigation above has clearly shown that, whatever the 
bath-mode spectrum is, no QPT would happen if the lo- 
cal field is absent in the SBM. 

In summary, we have presented analytical treatments 
to show reliably the impossibility of the QPT in the SBM 
without the local field. Since our treatments are inde- 
pendent of the bath type, the degrees of freedom of the 
bath and the calculation precision, we argue that the pre- 
vious conclusions drawn for the QPT happening in the 
Ohmic and sub-Ohmic SBM need more serious reexam- 
ination. The parity discovered in the present work is 
strongly relevant to the configuration of the model and 
the parity breaking can significantly change the dynamics 
of the model. In addition, the SBM has wide application 
ranging from the electron transfer in biomolecules [24J to 
the entanglement in quantum information science [l|, [2J . 
Therefor, we believe that our study is of general interest 



we first check the case (1). Since we require 
e - 2 £ fc 9 fc ^ o, the case J2{ n } c f n } D {o,n} = implies 
E{«} c { n } IL [( 2 9/c)" fc /v / ^] = 0, which is not satisfied 
unless cr i = 0. So the case (1) makes no sense and can 
be dropped. 

For the case (2), if we have Yl{ n } c {„} D {0,n}/ c {0} = 
- E{„} c { n } D {o,n}/ c {o} j£ 0, which leads to E+ 0} = E^ 0} , 
then we reach 
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Yj C {n} D {0,n} = — Yl C {n} D {0-,n}, 
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where Ylfn\=to * s ^ ne summation excluding the vacuum 
state of each bosonic mode. Since both sides of the above 
equation are of multivariate polynomials of q k , we may 
compare the terms of q k to determine if such an equation 
is satisfied. To this end, we further deduce the equation 
as 



^ e v,n 
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It is evident that the equation could not be satisfied since 
no constant term exists in the right side of the equator. 
This is contradictory to the above presumption E7~ -, = 

-ErQi ■ So we have proven the impossibility of degeneracy 
in the case of A ^ 0. 
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